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ABSTRACT 

We report the results of our third survey for formaldehyde (H2CO) 6 cm maser 
emission in the Galaxy. Using the Very Large Array, we detected two new H 2 CO 
maser sources (G23.01— 0.41 and G25.83— 0.18), thus increasing the sample of 
known H 2 CO maser regions in the Galaxy to seven. We review the characteristics 
of the G23.01— 0.41 and G25.83— 0.18 star forming regions. The H2CO masers 
in G23.01— 0.41 and G25.83— 0.18 share several properties with the other known 
H 2 CO masers, in particular, emission from rich maser environments and close 
proximity to very young massive stellar objects. 
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1. Introduction 

The first formaldehyde (H2CO) 6 cm maser was discovered in 1974 toward NGC7538 
(Downes & Wilson 1974; Forster et al. 1980); in the following ~ 28 years, H2CO masers 
were detected only toward two other regions in the Galaxy: Sgr B2 (Whiteoak & Gardner 
1983) and G29.96— 0.02 (Pratap et al. 1994). This low number of H 2 CO maser regions was 
unexpected (e.g., Forster et al. 1985; Gardner et al. 1986) given the detection of many H 2 CO 
maser spots in a single source (Sgr B2, Whiteoak & Gardner 1983; Mehringer et al. 1994), 
the widespread distribution of formaldehyde molecules as exemplified by Galactic H2CO 6 cm 
absorption (e.g., Watson et al. 2003; Araya et al. 2002; Downes et al. 1980; Dieter 1973), 
and the apparently common astrophysical conditions needed for the maser excitation if the 
masers were pumped by radio continuum radiation (see Pratap et al. 1992). 

The idea that H2CO masers may be pumped by background radio continuum was ini- 
tially proposed by Boland & de Jong (1981) to explain the maser in NGC7538 IRS1. How- 
ever, the low detection rate of new H 2 CO masers in dedicated surveys (Mehringer et al. 1995; 
Forster et al. 1985), the non-detection of maser emission from the H 2 CO 2 cm transition (e.g., 
Hoffman et al. 2003) and low emission measure (or non-detection) of radio continuum sources 
near several of the known H 2 CO masers, indicate that the pumping mechanism in most cases 
cannot be due to radio continuum excitation (e.g., see Araya et al. 2007c for the case of 
IRAS 18566+0408). The low detection rate of H 2 CO masers led several authors to speculate 
that formaldehyde masers are rare because specific and/or short-lived physical conditions 
may be needed for the excitation (e.g., Forster et al. 1985; Mehringer et al. 1995; Araya et 
al. 2007a). 

In an effort to understand the H 2 CO maser phenomenon and its place during the forma- 
tion of massive stars, we conducted two surveys for H 2 CO masers using Arecibo, the Green 
Bank Telescope, and the VLA0 that resulted in the detection of two new Galactic H 2 CO 
maser regions (IRAS 18566+0408, with emission first detected by Araya et al. 2004 and 
confirmed to be a maser by Araya et al 2005; and G23.71-0.20 with details in Araya et al. 
2006, from a survey later reported by Araya et al. 2007b). In this article we present the 
results of our third survey for H 2 CO masers. 



1 The Very Large Array (VLA) is operated by the National Radio Astronomy Observatory (NRAO), a 
facility of the National Science Foundation operated under cooperative agreement by Associated Universities, 
Inc. 
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2. Observations 
2.1. VLA Survey 

The initial observations were conducted in 2006 June with the VLA during the BnA 
— > B reconfiguration (9 syn ~ 1.5"). The details of the observations are reported in Table 1. 
A total of 14 pointing positions were observed toward massive star forming regions. In 
some cases, there were multiple radio continuum objects within the primary beam of a given 
position. Nine of the targets were selected based on their GBT and Arecibo spectra (Araya 
et al. 2002; Watson et al. 2003; Sewilo et al. 2004; Araya et al. 2007b) that showed 
profiles suggesting H2CO emission blended with absorption. In addition, given that H2CO 
6 cm masers appear to originate mostly from regions that also harbor H 2 and Class II 
CH3OH masers (e.g., Araya et al. 2005), we complemented our sample with 5 targets from 
the methanol and water maser catalog of Szymczak et al. (2005). The targets from the 
Szymczak et al. (2005) sample were observed ~40 min on-source. In the case of the targets 
that were selected based on single dish H 2 CO spectra, the time on-source ranged between 
10 to 100 min depending on the expected intensity of the H 2 CO emission line candidates. 
In Table 2 we list the observed targets, phase tracking center, center bandpass velocity, 
secondary (phase) calibrator used, and single-channel rms of the final data cubes. The data 
reduction and imaging were conducted using the NRAO package AIPS following the standard 
procedure for spectral line observations. 



2.2. High Spectral Resolution Observations 

Follow-up H 2 CO 6 cm observations of G23.01— 0.41 and G25.83— 0.18 were conducted 
on 2007 February 16 with the VLA in the D configuration. The goals of the observations 
were to confirm the H 2 CO 6cm emission detected in our survey (§3), and to study the 
H 2 CO maser line profiles with higher spectral resolution (1.83kHz, 0.11km s _1 ). The VLA 
correlator was used in the 1 IF mode with a bandwidth of 0.78 MHz (48 km s _1 ) and 512 
channels. The central bandpass velocity was set to 70 and 87km s^ 1 for the G23.01— 0.41 
and G25.83— 0.18 observations, respectively. The phase tracking center of the G23.01— 0.41 
observations was RA = 18 h 34 m 40!30, Decl. = -09°00'00"0 (J2000); the phase tracking center 
of the G25.83-0.18 observations was RA = 18 h 39 m 03^60, Decl. = -06°24'50'.'0 (J2000). The 
quasars 3C286 and J1822— 096 were used as primary and secondary calibrators, respectively. 
We assumed a flux density of 7.52 Jy for 3C286, and measured a flux density of 2.29 Jy for 
J1822— 096. In addition, 3C48 was observed to check the flux density calibration. Using 
3C286 as primary calibrator to bootstrap the flux density of 3C48, we measured a flux 
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density of 5.60 Jy for 3C48, which agrees within 3% with the expected value of 5.47 Jy. The 
calibration and imaging were conducted using the NRAO package AIPS. 

3. Results 

Of the 14 pointing positions observed with the VLA (and a total of 18 targets, see Ta- 
ble 2) we detected H 2 CO 6 cm emission toward two sources: G23.01— 0.41 and G25.83— 0.18 
(Figure 1). In both cases the emission was spatially unresolved (deconvolved size < 0.5") 
implying brightness temperatures greater than 10 4 K, which indicates maser emission. Thus, 
G23.01— 0.41 and G25.83— 0.18 are the sixth and seventh regions in the Galaxy where H 2 CO 
6 cm masers have been detected. 

Figure 1 shows the low velocity resolution spectra and peak channel images of the 
two new maser regions. The G23.01— 0.41 line is spectrally unresolved whereas the H 2 CO 
emission in G25.83— 0.18 has a double peaked profile. Figure 2 shows the line profiles of the 
high spectral resolution (VLA-D, §2.2) observations. H2CO emission was confirmed in both 
cases. In Table 3 we list the line parameters of the H 2 CO 6 cm masers from both observing 
periods. 

We detected radio continuum emission as well as H 2 CO 6 cm absorption toward several 
of the targets in the sample; a discussion of these data will be the topic of a future paper. 
Here we focus on the H 2 CO 6 cm maser emission, and discuss the H 2 CO 6 cm absorption 
only in the case of G25.83— 0.18 since the H 2 CO maser originates at the (projected) spatial 
center of a molecular clump traced by H 2 CO absorption (§4.2). 

4. Discussion 
4.1. G23.01-0.41 

4-1.1. H 2 CO 6cm Maser Emission 

H 2 CO 6 cm maser emission was detected toward G23.01— 0.41 with the VLA in the BnA 
— > B reconfiguration (Figure 1), and confirmed with VLA-D observations conducted approx- 
imately 8 months later (§2.2, Figure 2). The peak position, LSR velocity, and line width 
measurements are consistent in both runs (Table 3). If the VLA-D spectrum is smoothed 
to the channel width of the lower spectral resolution observations (VLA BnA^B, Table 3), 
the peak intensities are consistent within 2a, revealing no apparent variability of the line in 
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a ~ 8 month period. 

As mentioned in §3, the H 2 CO 6 cm maser in G23.01— 0.41 was not spectrally resolved 
in the VLA BnA^B observations (line width < 1.0km s _1 , Figure 1). In the high spectral 
resolution observations (Figure 2) the line was barely resolved, and appears to be the super- 
position of two components. The separation of the two possible overlapping lines is <0.4km 



4-1.2. Infrared Environment 

Far-IR emission is found at the location of the H2CO maser as shown by IRAS dataB 
However, the far-IR emission at the H2CO maser location is overlapped with several nearby 
sources. MSX 21.4/im data toward G23.01-0.41 is also affected by confusion due to blended 
emission with a source ~ 20" east of the H 2 CO maser position. At 70/xm, MIPSGAL 
Spitzer data show that the maser is coincident with a compact (FWHM < 30") and strong 
(saturated) FIR source. Since the emission is saturated, a reliable measurement of the 70 /zm 
size and position cannot be obtained; nevertheless, the peak position appears to be within 
~ 3" of the H2CO maser, i.e., certainly within the telescope angular resolution at 70 fim 
(~ 20"; Rieke et al. 2004). At the 24 /mi Spitzer MIPS band, the IR source coincident with 
the H 2 CO maser is substantially less affected by saturation and the center position can be 
more reliably determined; the 24 /zm source and the H2CO maser are less that 0.3" apart. 

The IR source coincident with the H2CO maser is detected in mid/near-IR observations 
from the Spitzer/IRAC GLIMPSE survey (Benjamin et al. 2003). Figure 3 (upper panel) 
shows a three color image (3.6/im blue, 4.5/xm green, 8.0/zm red) of the G23.01— 0.41 region. 
The Spitzer source at the H2CO maser location shows strong 4.5/xm excess indicative of 
shocked gas (e.g., Smith et al. 2006). There is also extended 4.5/mi excess emission west of 
the H 2 CO maser position (Figure 3, inset upper panel). 

Testi et al. (1998) detected a near infrared (NIR) source toward G23.01— 0.41, that 
was also detected by 2MASS. The 2MASS source is ~ 4" displaced from the H2CO maser 
position. Given this large angular offset, it is unlikely that the NIR source is responsible for 
the excitation of the masers in G23.01— 0.41 (§4.1.4). 



2 IRSA-IRAS Sky Survey Atlas, 



http://irsa.ipac.caltech.edu/Missions/iras.html 
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4-1.3. Radio Continuum and Thermal Molecular Lines 

We detected no 6 cm radio continuum toward the maser position to a 5a level of 
2.4 mJy beam -1 (6 syn = 1'.'6 x l'.'l). Using the VLA in the C configuration (6 sys ~ 1"), Codella 
et al. (1997) report an upper limit of the 1.3 cm radio continuum of 0.4mJybeam _1 (3a). 

Based on our VLA-D observations, the nearest 6 cm radio continuum source is approxi- 
mately 50" north of the maser position. The radio continuum source (a compact H II region 
that will be discussed in a future paper) has associated H 2 CO absorption at approximately 
the same velocity as the H 2 CO maser, indicating that the H 2 CO maser belongs to an ex- 
tended molecular cloud that shows clear evidence for the presence of more evolved massive 
stars. The LSR velocity of the H2CO absorption (75.0 km s _1 ) implies two possible kine- 
matic distances: 4.8 or 10.8 kpc. Codella et al. (1997) as well as Caswell & Haynes (1983) 
preferred the far kinematic distanced; however the distance ambiguity has not been resolved 
(Forster & Caswell 1999). For example, Harju et al. (1998) list the near kinematic distance 
for the source (see also Scoville et al. 1987, Pestalozzi et al. 2005). 

A number of thermal molecular lines have been detected toward G23.01— 0.41. Caswell 
et al. (2000) detected (quasi-)thermal emission of the CH 3 OH lines at 107.0 GH^| and 
156.6 GHz (see also Slysh et al. 1999). Harju et al. (1998) detected broad (line width > 
20 km s _1 ) SiO J=2-l and J=3-2 emission with the SEST telescope. CO and 13 CO emission 
was also detected in single dish surveys (Scoville et al. 1987; Jackson et al. 2006). 

Interferometric (Nobeyama Millimeter Array and Plateau de Bure Interferometer) ob- 
servations of 12 CO, 13 CO, and C 18 have been recently reported by Furuya et al. (2008) 
with an angular resolution ~6". Their 12 CO spectrum shows high velocity gas, in particular 
prominent red-wing emission that traces a massive (> 50 M©) molecular outflow, possibly 
located close to the plane of the sky. 

High density molecular gas in the region is evident from detection of CS (J = 1 — 0) 
and NH 3 (1,1) lines with the MIT Haystack 37 m telescope (Anglada et al. 1996). Single 
dish NH 3 (1,1), (2,2), and (3,3) observations were also conducted by Codella et al. (1997) 
with the Medicina 32 m telescope. Codella et al. (1997) derived a rotational temperature of 
13 K. 



3 Caswell & Haynes (1983; see also Caswell et al. 1995b, Anglada et al. 1996, Testi et al. 1998, Forster 
& Caswell 1999) report a far kinematic distance of 12.8 kpc instead of our value of 10.8 kpc because they 
assumed a Sun - Galactic Center distance of 10.0 kpc, whereas we assume 8.5 kpc (Brand & Blitz 1993). 

4 A maser component was also found superimposed on the thermal profile of the CH3OH 107.0 GHz line. 
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High angular resolution NH 3 observations were conducted by Codella et al. (1997) 
with the VLA (9 syn ~ 1"). The NH 3 (3,3) data show a compact (< 10") molecular core 
that has a SE-NW elongation (Figure 3, lower panel). Codella et al. (1997) derived the 
following parameters for the ammonia core: deconvolved angular diameter = 3'/3, T& = 58 K, 
M vir = 886 M©, and hh 2 = 6.9 x 10 6 cm~ 3 . The H 2 CO maser is coincident with the NH 3 
(3,3) peak emission (Figure 3, lower panel). 

Furuya et al. (2008) recently conducted a detailed multi-wavelength study of G23.01— 0.41 
at high angular resolution (< 10") using CO (see above), HNCO, and CH 3 CN. Their ob- 
servations confirm the presence of a hot molecular core characterized by a CH 3 CN rotation 
temperature of ~ 120 K and a core mass of ~ 380 M . The velocity distribution of the 
CH 3 CN gas is consistent with rotation of a molecular core oriented almost perpendicular to 
the outflow direction. The elongation and velocity gradient of the CH 3 CN core is parallel to 
the elongation of the NH 3 emission shown in Figure 3. 

4- 1.4- Other Astrophysical Masers 

A variety of masers have been detected toward G23.01— 0.41 in single dish surveys. 
Methanol masers have been found at 6.7 GHz (Menten 1991; Caswell et al. 1995b; Szymczak 
et al. 2002; see also catalogs by Pestalozzi et al. 2005 and Xu et al. 2003), 12 GHz (MacLeod 
et al. 1993; Caswell et al. 1995a; Blaszkiewicz & Kus 2004), 107.0 GHz (Caswell et al. 
2000), 44 GHz (Slysh et al. 1994), and 95 GHz (Val'tts et al. 2000). H 2 22 GHz and 
OH ( 2 n 3/2 J = 3/2 ground state) masers have also been reported (Szymczak et al. 2005; 
Szymczak & Gerard 2004; Caswell & Haynes 1983). 

High angular resolution observations of the CH 3 OH 6.7 GHz masers have been conducted 
(Caswell, unpublished data; see Caswell et al. 2000); the position of the CH 3 OH maser is 
coincident with that of the H 2 CO maser within 1" rms. Forster & Caswell (1989, see also 
Forster & Caswell 1999) conducted VLA observations of H 2 22 GHz (6 syn = 3"9 x 1"8) and 
OH 1665 MHz {6 syn = 6"5 x 1"0) masers in G23.01-0.41. The H 2 CO maser is located within 
2" of the OH and H 2 maser positions (see Figure 3 lower panel). 

4.1.5. The Nature of the G23.01-0.41 Star Forming Region 

The abundant available multi-wavelength data shows that G23.01— 0.41 is an active 
site of massive star formation. The H 2 CO maser is located at the peak of a molecular 
core traced by NH 3 (3,3) and CH 3 CN. The hot molecular core also harbors a variety of 
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molecular masers. Given the highly confused far IR environment, it is not possible to reliably 
measure the bolometric luminosity of the massive stellar object pinpointed by the H 2 CO 
maser (Figure 3); however, based on the available mid and far IR data, the upper limit 
of the luminosity is ~ 10 6 L Q . As in the case of the infrared source associated with the 
H 2 CO maser in IRAS 18566 + 0408 (Araya et al. 2007c), the H 2 CO maser in G23.01-0.41 is 
found toward a source with 4.5/im infrared excess, which likely indicates shocked gas in an 
outflow. Considering the presence of a hot molecular core at the location of the H 2 CO maser, 
the detection of several other maser species, evidence for outflow and shocked gas based on 
CO, SiO and 4.5/im infrared excess emission, absence of radio continuum, and evidence for 
rotation of the molecular core, the H 2 CO maser appears to pinpoint the location of a very 
young massive stellar object in a evolutionary stage prior to the formation of a radio bright 
ultra-compact H II region. The central object may still be undergoing accretion; further 
observations should be made to clarify this point. 



4.2. G25.83-0.18 

4-2.1. H 2 CO 6 cm Maser Emission 

An the case of G23.01— 0.41 (§4.1), maser emission in G25.83— 0.18 was first detected 
in the VLA survey (§2.1) and then confirmed with higher spectral resolution observations 
(VLA-D, §2.2) ~ 8 months later. The H 2 CO maser shows a double peaked line profile 
(Figures 1 and 2). The separation between the two peak components is 1.5km s -1 . No 
significant difference in the maser flux density was found between the two epochs. 

As reported in Table 3, the line profile of the high spectral resolution observations is well 
fit by the superposition of three Gaussian profiles; two broad components (0.7 and 0.9 km 
s _1 FWHM) and a narrow component (FWHM = 0.29km s _1 ). Whether the line profile is 
composed of only two non-Gaussian maser lines or three (or more) components is unclear. 
For example, the profile could be due to an asymmetric line (the blue-shifted component) 
overlapped with a Gaussian line (the red-shifted component, see Figure 2). 



4-2.2. Infrared Environment 

The H 2 CO maser is at the center of an infrared dark cloud as revealed by 8.0/xm 
Spitzer/GLIMPSE observations (see Figure 4). The H 2 CO maser is located between a source 
with strong 4.5/zm excess emission (green in Figure 4, upper panel inset) and a source brighter 
at 8/xm (red in Figure 4). The peak of the 4.5 fim source is offset (~ 3") from the location of 
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the H2CO maser. No infrared emission was detected with MSX toward the position of the 
H 2 CO maser; no 2MASS source was found coincident with the H 2 CO maser either. 

Data from Spitzer MIPSGAL at 70/xm reveal a strong (> 100 Jy) far IR source whose 
peak is within from the position of the H 2 CO maser. At low level, the 70/im emission 
is extended and comprises a neighboring infrared source; however, the core emission is com- 
pact (FWHM < 25"), i.e., close to the theoretical telescope diffraction limit. The IR source 
coincident with the H 2 CO maser is also detected in the 24/im MIPS band. 

4-2.3. Radio Continuum and Thermal Molecular Lines 

We detected no radio continuum toward the position of the H 2 CO maser to a level of 
2.8 mJy beam -1 (5a) in the VLA BnA — > B observations (9 syn = 1"6 x l'.'l); the 5a upper 
limit set by the VLA-D data is 5mJybeam _1 (see Figure 4, upper panel). The nearest radio 
continuum source to the H 2 CO maser is ~2' to the west (G25.80— 0.16; Figure 4, upper panel) 
which was also detected in the continuum at 8.64 and 6.67 GHz by Walsh et al. (1998). 

The H 2 CO maser is coincident with a mm and sub-mm core. Walsh et al. (2003) 
detected compact 450/im and 850/im emission; the H 2 CO maser is located within 4" of 
the peak position of the sub-mm source (the JCMT beam is approximately 8" and 15" at 
450 fim and 850/im, respectively). Hill et al. (2005) conducted SEST/SIMBA observations 
of G25.83— 0.18 and detected 1.2 mm emission (S u = 5.4 Jy); the FWHM of the mm source 
is 60". Based on the mm detection, Hill et al. (2005) report a mass of 2.8 x 10 3 M or 
8.8 x 10 3 M Q depending on the kinematic distance (see below). The H 2 CO maser is coincident 
with the 1.2 mm peak within ~ 1". 

We detected H 2 CO 6 cm absorption in the region (Figure 4, middle and lower panels, 
VLA-D observations). The main H 2 CO absorption clump is coincident with the infrared dark 
cloud, and shows a shell-like brightness distribution (Figure 4, middle paneljf). The H 2 CO 
maser is located at the projected spatial center of the shell and near the radial velocity edge of 
the H 2 CO absorption line (Figure 4, lower panel). The shell brightness distribution could be 
due to: 1.) H 2 CO in gas phase is less abundant in the inner regions of the shell, for example, 
due to chemical gradients, smaller total molecular density, or depletion (see Young et al. 
2004 for the case of pre-protostellar cores); and/or 2.) the excitation conditions for H 2 CO 



5 The shell-like distribution of H2CO absorption (Figure 4, middle panel) is not due to overlapping H2CO 
maser emission and absorption at the central position; the measured velocity of the H2CO maser is not in 
the velocity range used to obtain the H 2 CO absorption image. 
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absorption are less favorable in the inner regions of the clump (e.g., Zhou et al. 1990 explained 
the detection of an H 2 CO shell-like structure in B335 as a consequence of higher molecular 
density at the core that quenches the anomalous absorption). Observations of other H 2 CO 
transitions as well as other mid and high density tracers are needed to fully investigate 
the nature of the shell-like H2CO absorption source. Nevertheless, the recent high angular 
resolution detection of NH 3 at the position of the H 2 CO depression (Longmore et al. 2007, 
see below) suggests that the H 2 CO shell-like structure is due to high density (> 10 5 cm~ 3 ) 
molecular gas in the center of the molecular core that quenches H 2 CO anomalous absorption. 

We also detected H 2 CO absorption close to the radio continuum source (compare Fig- 
ure 4 upper and middle panels). The similar velocity of the H 2 CO absorption gas associated 
with the continuum source and with the infrared dark cloud indicates that both are part of 
the same star forming complex. Assuming that the LSR velocity of the H 2 CO absorption 
line traces the systemic velocity of the cloud, the two possible kinematic distances to this 
massive star forming region are 5.6 and 9.7 kpc. Since the H 2 CO maser is associated with 
an infrared cloud seen in absorption against the mid-IR galactic background (Figure 4), it 
is likely that the region is located at the near kinematic position. 

G25.83-0.18 has been detected in CH 3 CN, HCO+, and H 13 CO + with the Mopra tele- 
scope by Purcell et al. (2006). Based on the CH3CN data, they derived a rotation temper- 
ature of ~ 50 K, and found evidence for infalling motion based on HCO + data. In addition, 
Longmore et al. (2007) detected an optically thick NH 3 core at the position of the H 2 CO 
maser from ATCA observations. The line width of the NH 3 (1,1), (2,2), (4,4), and (5,5) 
transitions range between 6 and 30 km s _1 , indicating the possible presence of a molecular 
outflow partially traced by NH 3 . 

4-2.4- Other Astrophysical Masers 

Walsh et al. (1998) conducted ATCA observations of 6.7 GHz CH 3 OH masers in 
G25.83— 0.18; five maser components were found with LSR velocities between of 90.7 and 
99.0 km s _1 . The H 2 CO maser is located ~ 2" south of the CH 3 OH maser clump^l (Figure 4, 
upper panel, inset). 

Single dish observations of the Class II CH 3 OH 6.7 and 12 GHz masers were conducted 
by Blaszkiewicz & Kus (2004) with the Toruh telescope. The peak LSR velocity of the 6.7 
and 12 GHz masers were 91.3 and 90.7km s _1 , i.e., coincident in velocity with the H 2 CO 



The absolute position accuracy of the CH3OH masers is ~ 1", Walsh et al. (1998). 
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maser. 

Ellingsen (2005) conducted observations with the Mopra Telescope of the 95.1 GHz 
Class I CH 3 OH maser transition, and detected several maser lines (peak maser emission at 
90.2 km s _1 ) overlapped with a broad component. The broad line (4.3 km s -1 line width) 
could be due to thermal emission given that its peak velocity (94.2 km s _1 ) coincides with 
the H 2 CO absorption peak velocity (see §4.2.3). 

Other maser species have also been detected in single dish surveys. Szymczak & Gerard 
(2004) detected a OH 1667 MHz maser line with a peak LSR velocity of 92.9 km s^ 1 ; their 
figure A.l also shows detection of a possible high velocity (~ 120 km s -1 ) OH 1612 MHz 
line. Szymczak et al. (2005) conducted H2O 22 GHz observations with the 100 m Effelsberg 
telescope and detected several H 2 maser features within a velocity range of ~ 50 km s -1 , 
centered at 94.7 km s _1 . 

4.2.5. The Nature of the G25.83-0.18 Star Forming Region 

The coincidence of the H2CO maser with other maser species (in particular with Class 
II CH3OH masers), with an infrared dark cloud and molecular core, the absence of radio 
continuum emission, the detection of a mm, sub-mm and far-IR source, and the presence of 
4.5/xm excess emission toward the center of the infrared dark cloud, imply that G25.83— 0.18 
is a very young region of massive star formation in an evolutionary stage prior to the ultra- 
compact H II region phase. The precise location and luminosity of the protostar responsible 
for the excitation of the different maser species is unclear; higher sensitivity continuum 
observations are needed to reveal the position of the exciting source. G25.83— 0.18, as well 
as G23.01— 0.41, may be classified as Group 2 cores following the Longmore et al. (2007) 
nomenclature, i.e., warm NH 3 cores associated with CH 3 OH masers but no detectable radio 
continuum. 



4.3. G23.01-0.41 and G25. 83-0.18 with respect to the other known H 2 CO 

Maser Regions 

Including the two new masers reported in this work, H 2 CO 6 cm masers have been 
detected toward seven regions in the Galaxy, and in a total of 15 maser 'spots' (at 1" res- 
olution). The H 2 CO masers in G23.01— 0.41 and G25.83— 0.18 share similar characteristics 
with most of the other known H 2 CO maser regions, in particular: (1.) H 2 CO masers are 
found in regions that harbor a variety of other molecular masers (e.g., Hoffman et al. 2003; 
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Mehringer et al. 1994), (2.) the velocity difference of the H 2 CO masers with respect to the 
systemic velocity of the clouds is typically less than 6 km s _1 (e.g., Pratap et al. 1994; Araya 
et al. 2004), which suggests that H 2 CO masers do not originate in high velocity outflows, 
(3.) the flux density of the known masers is less than ~2 Jy (e.g., Hoffman et al. 2007), 
(4.) all known H2CO masers have been detected in regions of massive star formation (Araya 
et al. in prep.), (5.) excluding Sgr B2, most of the H 2 CO masers show double peaked 
profiles with separations smaller than 3 km s" 1 (e.g., NGC 7538 IRS1, Forster et al. 1985; 
IRAS 18566+0408, Araya et al. 2007d; see also review by Araya et al. 2007a), (6.) even 
though bright radio continuum sources may be found in the same star forming complexes, 
most H 2 CO masers are located toward sources characterized by weak (or no) compact radio 
continuum emission, typically undetected at a few mJy sensitivity levels (e.g., Araya et al. 
2005; this work), (7) excluding some of the masers in Sgr B2, the known H 2 CO masers 
appear to be associated with very young massive stellar objects (in an evolutionary phase 
prior to the formation of radio-bright ultra-compact H II regions) that have strong far IR 
emission and molecular core counterparts (e.g., §4.1 and 4.2). 

Many H 2 CO maser regions show evidence of outflows/jets and shocked gas based on 
4.5/zm excess emission from Spitzer/IRAC data (see, for example, Araya et al. 2007c in the 
case of IRAS 18566+0408). Further evidence comes from molecular data such as SiO, H 2 , 
and H 2 S (e.g., Zhang et al. 2007; Beuther et al. 2007b; Kraus et al. 2006; Gibb et al. 2004; 
Maxia et al. 2001; Harju et al. 1998), and radio continuum observations (e.g., Araya et al. 
2007c). However, the H 2 CO masers may not be directly associated with the shocked material 
(e.g., note the offset between the 4.5/xm excess source and the H 2 CO maser in G25.83— 0.18; 
Figure 4, upper panel). 

5. Summary 

We report the results of our third survey for H 2 CO 6 cm masers toward massive star 
forming regions in the Galaxy. The observations were conducted with the VLA toward 14 
pointing positions, and resulted in the detection of two new H 2 CO maser regions: G23.01— 0.41 
and G25.83— 0.18. Including the new detections, H 2 CO masers have been found toward a 
total of seven star forming regions in the Galaxy, four of them detected in our series of 
surveys (Araya et al. 2004, 2007b, and this work). 

The H 2 CO maser in G23.01— 0.41 is coincident with the center of a hot molecular core 
that shows evidence of a rotating torus perpendicular to a molecular outflow (Furuya et 
al. 2008). Excess in the 4.5/zm band from GLIMPSE Spitzer/IRAC observations reveals 
shocked gas in the region, possibly tracing the outflow. Active massive star formation is also 
evident from the detection of a number of molecular maser lines, including Class II CH 3 OH 
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masers. 

In the case of G25.83— 0.18, the maser is located at the center of an infrared dark cloud 
that harbors a clump of CH 3 OH 6.7 GHz masers, a high density molecular core traced by 
NH3 and CH3CN (Purcell et al. 2006), a mm, sub-mm and far IR counterpart. In addition 
there is evidence for molecular infall and outflow motions, as well as shocked gas traced by 
4.5 jum excess from Spitzer GLIMPSE data. 

In both cases the H 2 CO masers are found toward objects that have no compact 6 cm 
radio continuum emission at the few mJy/beam level. Given the evidence of active massive 
star formation in G23.01— 0.41 and G25.83— 0.18, the absence of radio continuum suggests 
that the regions are in a very early phase of massive star formation, prior to the develop- 
ment of radio-bright ultra-compact HII regions. H 2 CO 6 cm masers appear to preferentially 
pinpoint very young massive stellar objects that may be categorized as HMPOs (High-Mass 
Protostellar objects, Beuther et al. 2007a) or, in the Zinnecker & Yorke (2007) nomenclature, 
somewhere in the HDMC (Hot Dense Massive Core) or DAMS (disk-accreting main-sequence 
star) phase. 

E. A. is supported by a NRAO predoctoral fellowship. H. L. was supported by a 
postdoctoral stipend from the German Max Planck Society. L. O. was supported in part by 
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Table I. VLA H 2 CO Survey 



Parameter 


Value 


Date 


2006 Jun. 13 - 15 


Configuration" 


BnA^B 


vo (GHz) 6 


4829.6594 


IF Mode 


2IF (AD) 


BW (MHz) c 


1.56 


(km s- 1 ) 


97.0 


Au d (kHz) 


7.3 


(km s- 1 ) 


0.46 


Flux Density Calib. 


3C286 


Assumed S„ (Jy) 


7.52 


Phase Calib. 


J1832-105 


Measured Si, (Jy) 


1.26 


Phase Calib. 


J1824+107 


Measured Si, (Jy) 


0.80 


Phase Calib. 


J1950+081 


Measured Si, (Jy) 


1.07 



a Rcconfiguration from BnA to B array. 

b Weighted average frequency of the H2CO 
Jjf a K c = lio - 111 F = 2 - 2 and F = - 
1 hyperfine transitions given in Tucker et al. 
(1970). 

c Bandwidth per IF. 

d Spectral resolution. 
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Table 2. Observed Targets. 



Target 




a(2000)* 


<5(2000)* 




V LSR 


Phase Calib. 


rms' 






(h m s) 


(° 


' ") 




(km s- 1 ) 




(mjy) 


G13.13-0. 


.15 


18 


14 41.27 


-17 


37 


04. 


2 


40.0 


J1832-105 


7.6 


G13.21-0 


14 




t 




t 






40.0 


J1832-105 


17.4 


G21.57-0. 


03 


18 


30 36.50 


-10 


06 


44. 





110.0 


J1832-105 


5.7 


G22.05+0. 


.22 


18 


30 35.70 


-09 


34 


26. 





50.0 


J1832-105 


5.2 


G22.36+0. 


.06 


18 


31 44.84 


-09 


22 


08. 


5 


85.0 


J1832-105 


4.9 


G23.01-0. 


.41 


18 


34 40.60 


-09 


00 


28. 





75.0 


J1832-105 


6.2 


G22.97-0. 


.39 




* 




★ 






75.0 


J1832-105 


5.4 


G24.68-0. 


16 


18 


36 51.80 


-07 


24 


48. 


3 


50.0 


J1832-105 


4.6 


G25.79-0. 


.14 


18 


38 51.73 


-06 


25 


00 


6 


100.0 


J1832-105 


4.9 


G25.80-0. 


.16 




* 




* 






100.0 


J1832-105 


5.2 


G25.83-0. 


.18 




* 




* 






100.0 


J1832-105 


6.7 


G27.28+0. 


.15 


18 


40 33.90 


-04 


57 


18. 





35.0 


J1832-105 


5.6 


G30.87+0. 


.11 


18 


47 15.76 


-01 


47 


10. 


3 


90.0 


J1832-105 


10.1 


G33.92+0. 


11 


18 


52 50.27 


+00 


55 


29. 


6 


100.0 


J1824+107 


4.8 


G37.87-0. 


40 


19 


01 53.57 


+04 


12 


49. 


2 


60.0 


J1824+107 


5.6 


G39.10+0. 


.48 


19 


00 59.50 


+05 


42 


28. 





15.0 


J1824+107 


5.9 


G42.43-0. 


26 


19 


09 50.08 


+08 


19 


39. 


.7 


60.0 


J1950+081 


3.6 


G43.26-0. 


21 


19 


11 05.74 


+09 


05 


06. 





60.0 


J1950+081 


9.3 



Note. — (*' Phase tracking center. (°) Center bandpass LSR velocity. W Sin- 
gle channel rms of the final H 2 CO data cube, (t) G13.21-0.14 was included in the 
VLA primary beam of G13.13— 0.15. The reported rms was obtained toward the 
G13.21— 0.14 position after correcting for primary beam attenuation. M G22.97— 0.39 
was included in the VLA primary beam of G23.01— 0.41. The reported rms was ob- 
tained toward the G22.97— 0.39 position after correcting for primary beam attenua- 
tion. (*) G25.83-0.18 and G25.80-0.16 were included in the VLA primary beam of 
G25.79— 0.14. The reported rms values were obtained toward the G25.83— 0.18 and 
G25.80— 0.16 positions after correcting for primary beam attenuation. 
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Table 3. Line Parameters of H 2 CO 6 cm Masers 



Source 


a peafe (J2000) 
(h m s) 


<5peafc(J2000) l v ,peak 

(" ' ") (mjy/b) 


Vlsr 
(km s- 1 ) 


FWHM 

(km s- 1 ) 


VLA BnA — > B (Low Spectral Resolution) Observations 


G23.01- 
G25.83- 


-0.41* 
-0.18* 


18 34 40.294(0.004) 
18 39 03.627(0.003) 


-09 00 38.26(0.08) 48(6) 
-06 24 11.18(0.06) 85(7) 
45(7) 


73.5(0.4) 
90.2(0.4) 
91.4(0.1) 


<1.0 
<0.8 
1.8(0.3) 






VLA-D (High S 


Spectral Resolution) Observations 






G23.01- 
G25.83- 


-0.41* 
-0.18* 


18 34 40.29(0.06) 
18 39 03.69(0.07) 


-09 00 37.7(1.3) 53(9) 
-06 24 10.4(0.8) 102(9) 
52(9) 
44(9) 


73.60(0.09) 
90.21(0.02) 
91.64(0.08) 
90.6(0.1) 


0.4(0.2) 
0.29(0.04) 
0.9(0.2) 
0.7(0.3) 



Note. — Unless indicated otherwise, we list the intensity of the peak channel (the rms is listed as 
uncertainty), the LSR velocity of the peak channel (channel separation reported as uncertainty), and 
the FWHM (two times the channel separation is reported as uncertainty), (t) 9 syn = l'/Q x l'.'l, P.A. 
= —7°. * Osyn = 1"6 X l'.'l, P.A. = —18°. The line parameters of the weak component were derived 
from a Gaussian fit. (*) d sy n = 21"xl2", P.A. = -13°. W 9 syn = 23"xl5", P.A. = -8°. The line 
parameters were obtained from a three component Gaussian fit; the intensity and velocity of the two 
peak channels are: (113 mjy beam -1 , 90.2km s -1 ) and (66 mjy beam -1 , 91.7km s -1 ); see Figure 2. 



-09 00 35 
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37 
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Fig. 1. — Detection of H 2 CO 6 cm masers toward G23.01— 0.41 (rms = 6.2 mJy beam 1 ; 
e syn = 1% x l'/l, P. A. = -7°) and G25.83-0.18 (rms = 6.7 mJy beam" 1 ; 9 syn = 1% x l'/l, 
P. A. = —18°). The peak channel images are shown in contours and the spectra are shown 
in the insets. The observations were conducted in 2006 June (see §2.1), with a spectral 
resolution of 0.46 km s _1 . 
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Fig. 2. — High spectral resolution (0.11km s 1 ) VLA-D observations of the H 2 CO 6 cm 
masers in G23.01-0.41 and G25.83-0.18. The G25.83-0.18 H 2 CO spectrum has a clear 
double peaked profile; the line profile of the maser in G23.01— 0.41 also appears to be the 
superposition of 2 lines. The observations were conducted in 2007 February (§2.2). 
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18 34 40.4 40.2 40 



RA (J2000) 

Fig. 3.- Upper Panel: Spitzer/IRAC GLIMPSE three-color image of the G23.01-0.41 
region; the 3.6, 4.5, and 8.0/mi bands are shown in blue, green, and red, respectively. The 
location and FWHM upper limit of the 70 /mi MIPSGAL counterpart is shown with a ma- 
genta circle (see §4.1.2). A closeup of the G23.01— 0.41 region is shown in the inset. Note 
the strong 4.5 /mi excess toward the H 2 CO maser that extends westwards. The ellipse rep- 
resents the approximate size and position angle of the NH 3 core shown in the lower panel 
(size, position and orientation of the ellipse are from a fit of the 12 mJybeam -1 contour 
of the Codella et al. (1997) map). The location of the H2CO maser is explicitly marked 
in all panels. Lower Panel: Zoomed-in view of G23.01— 0.41. The contours show the NH 3 
(3,3) emission imaged with the VLA by Codella et al. (1997; the original VLA image was 
convolved with a Gaussian profile whose FWHM is shown in the upper-left corner of the 
figure). The location of H 2 (22 GHz) and OH (1665 MHz) masers from Forster & Caswell 
(1999) are marked with squares and diamonds, respectively. The absolute astrometry of the 
OH and H 2 maser positions is ~ C/5 rms (see Forster & Caswell 1999). 




H 2 CO Absorption 



8 39 08 06 04 02 00 38 58 

RIGHT ASCENSION (J2000) 
Peak Absorption = -594.6 mjy/b 
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Fig. 4. — Upper panel: We show in contours the radio continuum (6 cm) image of the 
G25.83-0.18 region (rms = lmJybeam -1 ; 6 syn = 23"0 x 15"0, P.A. = -8°) superimposed 
on a Spitzer/IRAC GLIMPSE three-color image of the region (3.6/im blue, 4.5yum green, and 
8.0/iHi red). The location and FWHM upper limit of the 70 /iin counterpart is shown with 
a magenta circle (see §4.2.2). A zoomed-in view of the H 2 CO maser region is shown in the 
inset; the position of CH 3 OH 6.7 GHz masers are also shown (Walsh et al. 1998; absolute 
position accuracy ~ 1"). Middle Panel: Sum of the channel maps in the 92.4 to 95.5km 
s _1 velocity range. The position of the H 2 CO maser is shown with a '©' symbol in the upper 
and middle panels; the orientation of the symbol represents the P.A. of the 9 syn and its size 
is five times the size of the 8 syn (§2.1, Table 3). In the inset of the upper panel, the size of the 
symbol is that of the VLA synthesized beam. Lower Panel: Spectrum of the G25.83— 0.18 
maser region obtained by integrating the brightness distribution in a 85 // x90 // area centered 
at the H 2 CO maser position. Note the maser emission at ~ 90 km s -1 . 



